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ABSTRACT: 

Bioenergetic  modeling  provided  insight  into  the  potential  for  fishes  to  impact  zebra  mussel  Dreissem  polymorpha 
populations  by  predation  across  nver  and  lake  systems  in  eastern  North  America  (13°  latitude  x  12°  longitude). 

ood  consumption  modeling  suggested  that  fishes  in  southern  latitudes  consumed  up  to  100  percent  more  food  than 
those  m  northern  systems.  Much  variation  in  food  consumption  resulted  from  differing  water  temperature  regimes 
which  presumably  dictated  tfie  bioenergetic  demands  of  fishes.  The  potential  for  fishes  to  impact  zebra  mussels  was 
also  mfluenced  by  the  standing  crop  of  different  species  and  the  fish  community  as  a  whole.  Fish  community 
structures  followed  predictable  patterns  along  both  longitudinal  and  latitudinal  gradients.  Multivariate  analyses 
m  icated  a  tendency  for  central  and  southern  U.S.  systems  to  contain  greater  standing  crops  or  biomasses  of  fishes 
likely  to  consume  zebra  mussels.  Though  analyses  indicated  northern  systems  had  greater  biomasses  of  some 
potential  zebra  mussel  predators,  net  biomass  increases  of  these  species  was  not  great  enough  to  offset  the  decreases 
m  food  consumption  because  of  cooler  annual  water  temperatures.  This  study  generally  supports  the  premise  that 
fishes  in  more  southern  (including  central)  U.S.  waters  have  inherently  greater  potential  to  impact  zebra  mussels 
because  of  comm^ity  composition  and  bioenergetics.  Our  simulations  provide  a  partial  explanation  of  why  zebra 

mussel  mvasions  have  not  been  as  rapid  and  widespread  in  southern  U.S.  waters  compared  to  the  Great  Lakes 
region. 
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1  Introduction 


Zebra  mussels  Dreissena  polymorpha  invaded  North  America  in  the  mid- 
1980s  (Nalepa  and  Schloesser  1993).  They  are  perhaps  the  most  influential 
animals  to  invade  North  American  fresh  waters  (Thorp,  Delong,  and  Casper 
1998),  having  effects  across  numerous  ecosystem  levels  (Richardson  and  Bartsch 
1997;  Wyer  et  al.  2000;  Idrisi  et  al.  2001).  Zebra  mussels  were  first  collected  in 
Lake’st.  Clair  in  1986  (Ram  and  McMahon  1996),  spread  rapidly  through  the 
Great  Lakes  by  the  early  1990s  (Mills  et  al.  1993),  and  then  southward  through 
the  Illinois,  Mississippi,  Ohio,  Cumberland,  Tennessee,  Arkansas,  and 
Atchafalaya  Rivers  by  the  mid-1990s  (Allen,  Thompson,  and  Ramcharan  1999). 

Zebra  mussels  occur  patchily  in  most  systems  with  the  greatest  densities 
occurring  on  harder,  exposed  substrates.  Densities  commonly  range  from  2,000 
to  100,000  mussels-m'^  (Nalepa  and  Schloesser  1993),  though  densities  up  to 
340,000  mussels  m'^  have  been  reported  from  western  Lake  Erie  (Leach  1993). 
Greater  densities  have  been  reported  anecdotally  from  other  systems.  At  the 
southern  end  of  their  range,  the  lower  Mississippi  River,  densities  are  merest, 
typically  on  the  order  of  hundreds  of  mussels  m'^  (Personal  Communication 
D.  Miller,  U.S.  Army  Corps  of  Engineers  (USACE)). 

Zebra  mussels  thrive  in  northern  United  States  (U.S.)  waters  where  water 
temperatures  usually  do  not  exceed  27  to  28  °C  (McMahon  1996).  Aldridge, 
Payne,  and  Miller  (1995)  found  high  survival  of  zebra  mussels  at  temperatures  up 
to  32  °C  in  laboratory  experiments.  Allen,  Thompson,  and  Ramcharan  (1999) 
inferred  that  zebra  mussel  reproduction  and  survival  would  likely  be  adversely 
affected  by  the  higher  water  temperatures  present  in  most  southern  U.S.  river 
systems.  However,  southern  zebra  mussel  populations  are  firmly  established  ^ 
despite  initial  predictions  that  summertime  temperatures  would  preclude  sitrvival 
(Allen,  Thompson,  and  Ramcharan  1999). 

Another  factor  inhibiting  zebra  mussel  expansion  in  southern  waters  may  be 
predation  by  native  fishes.  Eggleton  (2001)  reported  extensive  (up  to  26  percent 
of  annual  food  consumption)  use  of  zebra  mussels  by  blue  catfish  Ictalurus 
fiircatus  in  the  main  channel  of  the  lower  Mississippi  River.  Other  regional 
reports  indicate  blue  catfish  and  freshwater  drum  Aplodinotus  grunniens  in  the 
Arkansas  River,  Arkansas  (Magoulick  and  Lewis  2002);  common  carp  Cyprinus 
carpio,  other  catfishes  Ictalurus  spp.,  and  freshwater  drum  in  the  Cuniberland 
River,  Tennessee  (Personal  Communication,  D.  Hubbs,  Tennessee  Wildlife 
Resorlrces  Agency);  and  blue  catfish,  river  carpsucker  Carpiodes  carpio,  and 
smallmouth  buffalo  Ictiobus  bubalus  in  the  Ohio  River  (Thorp,  Delong,  and 
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Casper  1998)  also  consume  zebra  mussels.  Similar  reports  also  exist  from 
northern  U.S.  waters  (e.g.,  French  and  Bur  1993;  Tucker  et  al.  1996;  French 
1997;  Marsden  1997;  Morrison,  Lynch,  and  Dabrowski  1997),  principally  with 
common  carp  and  freshwater  drum.  Molloy  et  al.  (1997)  reported  that  up  to 
38  species  of  fish  eat  zebra  mussels,  though  some  species  may  be  less  adept 
because  of  gape  limitations  (e.g.,  French  and  Love  1995).  However,  they  added 
that  the  few  studies  conducted  suggest  that  fish  predation  does  not  significantly 
affect  zebra  mussel  populations. 

There  are  several  possible  reasons  for  lack  of  predatory  control  in  northern 
U.S.  waters.  Physical  conditions  in  many  northern  rivers  and  lakes  mimic  those 
in  Eurasia,  the  native  range  of  zebra  mussels  (Mackie  and  Schloesser  1996),  and 
are  conducive  to  their  survival  and  proliferation.  Many  northern  U.S.  waters  also 
have  lower  population  densities  of  large-bodied  omnivorous  fishes  (e.g., 
catostomids,  catfishes,  and  freshwater  drum)  that  could  consume  signifiknt 
amounts  of  zebra  mussels.  This  attribute,  coupled  with  cooler  water  temperatures 
and  lower  bioenergetic  demands  of  fishes  (e.g.,  growth  and  metabolism)  in  most 
northern  waters,  suggests  that  fishes  will  consume  less  food.  Thus,  fish 
communities  in  many  northern  U.S.  rivers  and  lakes  may  have  less  potential  to 
impact  zebra  mussels.  On  the  other  hand,  southern  waters  are  generally  warmer 
and  often  have  greater  standing  crops  of  potential  mussel  predators 
(molluscivores)  that  include  ictalurid  catfishes,  common  carp,  catostomid 
suckers,  and  freshwater  drum.  Growing  seasons  are  longer  in  southern  systems 
and  allow  resident  fishes  to  exhibit  high  growth  rates  and  attain  large  sizes. 

Thus,  fishes  in  southern  U.S.  waters  may  have  a  competitive  edge  in  controlling 
or  restricting  zebra  mussels  not  possible  in  northern  systems 

The  regulation  or  control  of  zebra  mussels  and  other  exotic  molluscs  by 
fishes  through  predation  has  not  been  widely  assessed.  French  and  Bur  (1993) 
provided  a  review  on  the  potential  of  several  North  American  fishes  to  utilize 
zebra  mussels  based  on  the  mouthpart  morphology.  Thorp,  Delong,  and  Casper 
(1998)  conducted  a  series  of  enclosure-exclosure  experiments  in  the  Upper 
Mississippi  and  Ohio  rivers.  They  found  fishes  had  a  greater  impact  when  zebra 
mussels  occurred  at  low  densities,  but  fish  predation  was  still  insufficient  to 
impact  mussels.  Conversely,  Robinson  and  Wellborn  (1988)  concluded  that 
fishes  in  a  Texas  reservoir  had  great  capacity  to  impact  the  invading  Asiatic  clam 
Corbiculafluminea,  exhibiting  as  much  as  a  29-fold  effect  on  clam  abundance  in 
cages  where  they  had  access.  Thus,  the  few  field  studies  conducted  suggest  the 
ability  of  fishes  to  impact  exotic  mussels  may  vary  among  systems  and  be  related 
to  latitude. 


Enclosure-exclosure  experiments  provide  insight  into  predatory  mechanisms 
of  many  animals.  However,  the  field  of  bioenergetics  offers  another  approach  to 
assess  the  potential  for  fishes  to  impact  prey  populations  (e.g.,  Stewart,  Kitchell, 
and  Crowder  198 1).  Simple  bioenergetic  techniques,  termed  food  consumption 
models,  estimate  the  daily  food  consumption  (ration)  rate  of  fishes  as  a  fimction 
of  the  amount  of  food  observed  in  the  stomach,  gastric  evacuation  rate,  and  water 
temperature  (e.g.,  Eggers  1977;  Elliott  and  Persson  1978;  Boisclair  and  Leggett 
1988).  Other  approaches  estimate  population-level  food  consumption  as  a 
function  of  fish  production  rates,  gross  conversion  efficiencies,  or  production- 
biomass  ratios  (e.g.,  Eck  and  Wells  1983;  Roell  and  Orth  1993).  These 
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approaches  are  termed  production-  or  hiomass-based  models.  More  elaborate 
techniques,  termed  mass-balance  models,  relate  the  feeding  and  growth  rates  of 
fishes  through  construction  of  energetic  mass  budgets  that  balance  energy 
acquisition  (i.e.,  food  consumption)  with  energy  expenditures  (i.e,,  metabolism, 
growth)  and  losses  (i.e.,  waste)  (Tytler  and  Calow  1985),  Recent  models 
developed  for  this  purpose  (e,g.,  the  ‘Wisconsin  model’  of  Kitchell,  Stewart,  and 
Weininger  1977)  use  a  series  of  sophisticated  species-specific  mass-balance 
algorithms  applied  to  individual  fishes  and  extrapolated  to  whole  populations 
(e.g,,  Kitchell,  Stewart,  and  Weininger  1977;  Lyons  and  Magnuson  1987;  Rabom 
2000).  The  mass-balance  approach  has  become  the  most  common  in  North 
American  fisheries  research  (Ney  1993),  but  it  is  restricted  to  species  for  which 
physiological  parameters  required  by  the  model  have  been  estimated.  Parameters 
available  (up  to  30  per  species)  are  almost  exclusively  limited  to  fishes  common 
to  the  Great  Lakes  region  (Hanson  et  al.  1997).  However,  all  three  types  of 
bioenergetic  models  could  be  used  to  estimate  potential  impacts  to  zebra  mussel 
populations  across  a  latitudinal  gradient. 

This  study  uses  a  bioenergetic-based  approach  to  investigate  latitudinal 
effects  on  food  consumption  by  fishes.  We  hypothesize  that  these  latitudinal 
relationships  coupled  with  regional  differences  in  fish  faunas  and  biomasses  may 
be  used  to  infer  the  relative  potential  for  fishes  to  impact  zebra  mussels.  Because 
physiological  data  needed  to  parameterize  a  mass-balance  model  were  largely 
imavailable,  our  evaluation  focused  on  the  food  consumption  and  a  biomass- 
based  models. 

Our  first  objective  was  to  simulate  the  annual  food  consumption  of  selected 
fishes  in  eight  representative  U.S.  river  and  lake  systems  infested  with  zebra 
mussels.  Daily  food  consumption  was  estimated  on  an  individual-fish  basis 
using  a  conventional  food  consumption  model  that  integrated  a  generalized 
gastric  evacuation  model  and  average  daily  water  temperature.  This  modeling 
yielded  species-specific  estimates  of  daily  food  consumption  for  each  day  of  a 
typical  year,  which  could  then  be  summed  to  determine  annual  food 
consumption.  The  second  objective  involved  simulating  population-level 
consumption  based  on  biomass  for  the  same  species  in  each  system.  A 
multivariate  spatial  analysis  of  fish  communities  across  geographical  regions  was 
also  conducted  to  augment  biomass  analyses.  To  accomplish  each  objective,  we 
had  to:  (a)  parameterize  representative  models  for  each  species  and  system  by 
review  of  pertinent  literature;  (b)  compare  predicted  food  consumption  among 
species  and  across  systems  that  represented  different  latitudesAongitudes  and 
thermal  regimes;  and  (c)  use  sensitivity  analyses  to  assess  the  robustness  of 
predictions  and  the  relative  importance  of  different  model  parameters.  This 
overall  approach  differs  from  previous  studies  (e.g.,  Robinson  and  Wellborn 
1988  and  Thorp,  Delong,  and  Casper  1998)  but  may  be  useful  to  resource 
managers  involved  with  forecasting  future  zebra  mussel  impacts  on  aquatic 
ecosystems. 
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2  Methods  and  Materials 


Study  Areas 

We  selected  eight  U.S.  river  and  lake  systems  that  were  infested  with  zebra 
mussels  for  which  adec^uate  thermal  data  existed.  Selected  study  systems  were 
located  along  a  13-degree  latitudinal  gradient  and  a  12-degree  longitudinal 
gradient  from  Louisiana  north  to  Minnesota  and  east  to  Lake  Erie  (Figure  1). 
These  systems  included  northern  Lake  Michigan,  Upper  Mississippi  River 
(UMR)  Pools  4,  13,  and  26;  Fox  River;  eastern  Lake  Erie;  the  upper  Ohio  River; 
and  the  Lower  Mississippi  River  (LMR)  (Table  1 ). 
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1  Table  1 

1  River  and  Lake  Systems  Selected  for  Food  Consump 
1  Simulations 

ition  Modeling 

System 

Location 

Latitude  fN) 

Longitude  ("W) 

Lower  Mississippi  River 

Mississippi-Arkansas 

32.3089 

90.908 

Upper  Miss  River  Pool  26 

lllinois-Missouri 

38.9276 

90.2803 

1  Ohio  River 

Ohio-West  Virginia 

39.1072 

81.7429 

Upper  Miss  River  Pool  13 

lllinois-lowa 

41.9594 

90.1565 

Lake  Erie 

New  York 

42.8613 

78.9215  1 

Fox  River 

Wisconsin 

44.3075 

88.2185  1 

Upper  Miss  River  Pool  4 

Minnesota-Wisconsin 

44.4202 

92.1183  1 

I  Lake  Michigan 

Wisconsin 

45.3300 

86.4200  I 

Selected  Fish  Species 

Molloy  et  al.  (1997)  reported  up  to  38  species  of  fish  consume  attached 
(nonveliger  or  adult)  zebra  mussels.  However,  geographic  variation  suggested 
that  our  best  candidates  were  large,  abundant,  and  geographically  widespread 
omnivorous  species.  Given  these  criteria,  we  judged  the  best  candidates  included: 
common  carp,  freshwater  drum,  channel  catfish  Ictalurus  punctatus,  blue  catfish, 
lake  sturgeon  Acipenser  julvescens,  spotted  sucker  Minytrema  melanops,  river 
carpsucker,  buffalofishes  Ictiobus  spp.,  and  redhorse  suckers  Moxostoma  spp.. 
Because  the  latter  two  groups  contained  several  species  of  similar  morphology 
and  overlapping  distributions,  collective  grouping  at  the  genus  level  was  done  in 
lieu  of  assessing  each  species  individually.  Eight  of  the  nine  species/genera  were 
modeled  using  the  biomass-based  approach — ^lake  sturgeon  were  excluded 
because  too  few  biomass  estimates  were  available  for  meaningful  analysis. 


Food  Consumption  Model 

Using  the  food  consumption  approach,  only  common  carp,  freshwater  drum, 
channel  catfish,  and  lake  sturgeon  were  modeled  because  of  the  extensive  time 
required  to  run  individual  simulations.  Given  that  only  the  growth  variable 
differs  among  species  in  the  food  consumption  model,  we  did  not  believe  running 
the  full  complement  of  simulations  (i.e.,  all  systems,  ages,  and  growth  scenarios) 
for  all  species  would  be  worthwhile. 

a.  Approximation  of  growth.  Estunates  of  fish  size  and  growth  were  needed 
to  simulate  annual  food  consumption.  Modeling  of  daily  food 
consumption  incorporated  estimates  of  stomach  fullness  expressed  as  a 
proportion  of  total  weight.  Because  we  needed  generalized  estimates  of 
fish  growth  across  a  large  geographic  area,  we  used  length-at-age  data 
reported  by  Carlander  (1969, 1977, 1997).  Data  sets  from  diverse 
geographic  areas  were  used  in  the  development  of  generalized  length-at- 
age  curves  that  were  then  converted  to  weight-at-age  curves  (Table  2) 
using  weight-length  relations  reported  in  Carlander  (1969, 1977, 1997). 
Mean  weight-at-age  curves  were  fiien  used  to  depict  generalized  annual 
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growth  curves  for  each  species.  In  addition,  given  that  fish  growth  can 
be  highly  system-specific  (DeVries  and  Frie  1996),  fast-growth  (mean 
weight-at-age  +  20  percent)  and  slow-growth  (mean  weight-at-age  - 
20  percent)  scenarios  were  also  modeled  to  provide  robustness. 
Simulations  were  conducted  for  three  age  intervals  for  each  species:  ages 
3  to  4,  5  to  6,  and  7  to  8. 


Table  2 

Generalized  Weight-Length  Regression  Equations  and  Welght-At-Age  Estimates  Used 
for  Each  Species  in  this  Evaluation 


Species 

Generalized  Weight- 
Length  Equations 

Estimated  Weights 

From  Ages  1-10 

Weight 

Interval 

Age  3-4 
(wt.  gain) 

Weight 

Interval 

Age  5-6 
(wt.  gain) 

Weight 

interval 

Age  7-8 
(wt.  gain) 

Common 

carp 

-4.489+2.832Log,o(TL) 

62,311,722, 1120, 1501, 
1918,  2172, 2624,3000, 
3423 

722-1120  (398) 

1501-1918 

(417) 

2172-2624 

(452) 

Freshwater 

drum 

-5.419+3.204Log,o(TL) 

29,  219, 492, 720,  848, 
1639, 1803, 2166,  2575, 
3035 

492-720  (268) 

848-1639  (791) 

1803-2366 

(563) 

Channel 

catfish 

-5.508+3.290Log,o(TL) 

13,  83,219,  439,  756, 

1226, 1653,  2324,  3904, 
4780 

219-439  (220) 

756-1226  (470) 

1653-2324 

(671) 

Lake 

sturgeon 

-3.482+2.837Log,o(FL) 

1816,4709,5511,9096, 
16670, 19747,  21095, 
23737,  28343,  40489, 
70609 

5511-9096 

(3585) 

16670-19747 

(3077) 

21095-23737 

(2642) 

Note:  Weight  data  were  used  to  generate  estimates  of  stomach  fuilness  for  each  species  during  simulations. 
TL~to^length  In  mm,  FL=fork  length  In  mm,  all  weights  in  grams  wet-weight. 


b.  Model  parameterization.  Our  goal  was  to  estimate,  mean  daily  food 
consumption  for  each  target  species  in  each  river  or  lake  system.  Each 
simulation  encompassed  a  365-day  period  beginning  on  July  1  (day  1) 
and  ending  on  June  30  (day  365).  Simulations  were  accomplished  using 
a  conventional  food  consumption  model,  a  model  for  generalized 
estimates  of  fish  gastric  evacuation  rate,  and  average  daily  water 
temperatures.  The  365  estimates  of  mean  daily  ration  were  then  summed 
to  determine  annual  food  consumption  for  a  typical  individual  of  each 
species  (i.e.,  common  carp,  freshwater  drum,  channel  catfish,  and  lake 
sturgeon)  in  each  system. 

Food  consumption  models  estimate  a  daily  (i.e.,  24-hr)  ration  by  integrating 
the  observed  amount  of  food  in  the  stomach  at  a  point  in  time  with  estimated 
instantaneous  rate  of  gastric  evacuation.  It  was  assumed  that  fishes  fed  randomly 
throughout  a  24-hr  cycle.  We  used  the  Eggers  food  consumption  model  (Eggers 
1977,  1979): 

Cj=2AxS,>^R  /j\ 
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where 


S,  =  weight-specific  amount  of  food  in  the  fish  stomaeh  determined  from 
stomach  fullness  (expressed  as  g  wet-weight-g  fish  *) 

R  =  instantaneous  evacuation  rate 
24  =  hours  in  a  day. 

Instantaneous  rates  of  gastric  evacuation  were  estimated  as  a  fimction  of 
water  temperature  and  meal  size  using  a  generalized  model  developed  for  fishes 
(He  and  Wurtsbaugh  1993): 

R  =  (2) 


where 


T  =  daily  mean  water  temperature  (°C) 

MS  =  estimated  meal  size  (g  wet-weight). 

Meal  size  was  approximated  based  on  a  mean  stomach  fullness  (i.e.,  weight 
of  stomach  contents  /  weight  of  the  fish)  of  1  percent  of  the  fish  total  weight  with 
a  possible  range  of  0  to  4.  These  ranges  of  stomach  fullness  were  based  on 
estimates  provided  in  the  literature  for  several  species  of  small  and  large-bodied 
fishes  (e.g.,  Armstrong  and  Brown  1983;  DeMauro  and  Miranda  1990;  Keast  and 
Fox  1992;  Driscoll  and  Miranda  1999).  Stomach  fullness  means  and  ranges  for 
two  of  the  target  species  collected  from  the  lower  Mississippi  River  (almost 
2,000  fishes  collected  over  2  years)  were  also  considered  (Eggleton,  M.  A. 
unpublished  data).  The  generalized  evacuation  model  of  He  and  Wurtsbaugh 
(1993)  was  based  on  121  published  gastrie  evacuation  estimates  that  included 
22  species  of  freshwater  fishes  (weights  10  to  1,300  g  wet-weight)  at  water 
temperatures  ranging  from  2  to  27  "C;  however,  none  of  our  target  species  were 
included  in  its  development.  Temperature  ranges  for  some  of  our  study  systems 
also  approximated  the  upper  and  lower  limits  of  temperatures  used  by  He  and 
Wurtsbaugh  in  their  model.  Thus,  to  add  robustness  to  our  daily  food 
consumption  predictions,  we  made  1,000  estimates  of  evacuation  rate  for  each 
day  of  each  simulation  using  random  estimates  of  the  model  coefficients,  and 
meal  size.  Random  estimates  of  the  evacuation  model  coefficients  were  selected 
within  triangular  eonstraints  of  the  90  percent  confidence  intervals  for  each 
coefficient  (Miranda,  Hargreaves,  and  Rabom  2001).  Similar  estimates  of  meal 
size  were  made  using  a  mean  of  1  percent  and  a  range  of  0  to  4  percent.  Because 
daily  water  temperatures  for  each  system  were  daily  means  from  at  least  several 
years  of  data,  estimates  were  believed  reasonably  precise  and,  therefore,  fixed  for 
each  day  during  simulations. 


Biomass-based  Model 

We  estimated  population-level  annual  food  consumption  for  eight  fish 
species  in  each  river  and  lake  system  as  a  function  of  production  rates  (derived 
from  biomass  data),  gross  conversion  efficiencies,  and  production-biomass  ratios. 
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Biomass  data  were  obtained  from  a  data  set  containing  estimates  for  all  fish 
species  from  396  U.S.  reservoirs.  Thirty-two  additional  biomass  estimates  were 
obtained  for  the  Upper  and  Lower  Mississippi  Rivers.  Reservoir  estimates  were 
compiled  primarily  from  state  fish  and  game  agencies  from  Florida  north  to 
Maryland  and  the  Ohio  River  and  west  to  Kansas  and  Texas.  River  biomass 
estimates  were  compiled  from  Federal  Government  Agencies,  namely  the 
U.S.  Army  Corps  of  Engineers  for  the  Lower  Mississippi  River  (Cobb  et  al. 
1984)  and  the  U.S.  Fish  and  Wildlife  Service  -  Upper  Mississippi  River 
Conservation  Commission  for  the  Upper  Mississippi  River  (Pitlo  1987). 

a.  Model  parameterization.  Using  these  data,  we  constructed  least-squares 
linear  regression  models  of  log,,  fish  biomass  versus  latitude,  longitude, 
and  latitude-longitude  interaction  for  each  target  species.  The  predicted 
biomass  of  each  study  system  was  used  to  estimate  the  expected  average 
biomass  such  that: 


expected  log^  (fish  biomass)  =  a  + 

b(latitude)+c(longitude)+d(interaction)  (3 

where  a,  b,  c,  and  d  were  estimated  parameters. 

The  product  of  the  expected  biomass  for  a  particular  system  and  the  range  of 
expected  production-biomass  (P-B)  ratios  for  warmwater  fishes  from  the 
literature  (mean  2.5,  range  0.5  to  4.5  (Neves  1981))  yielded  the  estimated  annual 
production  for  each  species  in  that  system.  Thus, 


P  =  BP-B 


(4) 


where 


B  =  predicted  biomass  in  kg-ha'’ 

P-B  =  production-biomass  ratio 

Dividing  the  predicted  annual  production  by  the  expected  range  of  gross 
conversion  efficiencies  (GCE)  for  fishes  (mean  0.1,  range  0.01  to  0.2;  Brett  and 
Groves  1979)  yielded  the  estimated  annual  food  consumption  for  each  species  in 
a  particular  system.  Thus, 

C,=P/{GCE)  (5 


where 


C,  -  population-level  annual  food  consumption  in  kg-ha'' -yf' 
P  =  annual  fish  production  in  kg-ha"' 

GCE  =  growth/consumption  (gross  conversion  efficiency) 
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Estimates  of  P-B  ratio  and  GCE  were  not  available  for  the  systems  or  species 
of  interest.  However,  individual  estimates  and  generalized  ranges  for  other 
species  and  systems  were  available  in  the  primary  literature.  Thus,  to  assure 
robustness  to  our  predictions,  5,000  random  estimates  of  predicted  biomass  were 
generated  for  each  species  and  system  within  triangular  constraints  of  the 
90-percent  confidence  intervals  of  the  log«(biomass)-latitude/longitude  regression 
models.  In  addition,  5,000  random  estimates  of  P-B  ratio  were  generated  using  a 
mean  of  2.5,  a  minimum  of  0.5,  and  a  maximum  of  4.5;  similarly,  5,000  random 
estimates  of  GCE  were  generated  using  a  minimum  of  1  percent,  mean  of 
10  percent,  and  a  maximum  of  20  percent.  The  product  of  each  set  of  biomass 
estimates,  P-B  ratios,  and  GCE  values  yielded  5,000  estimates  of  population- 
level  food  consumption  for  each  species  and  system.  Plotting  the  frequency  of 
these  expected  values  for  each  river  or  lake  system  simulated  a  distribution  of 
expected  or  potential  consumption  values  for  each  species  in  each  system  from 
which  means  and  standard  errors  could  be  estimated.  Given  the  randomization 
process,  this  approach  provided  a  reasonably  robust  yet  conservative  estimation 
of  potential  annual  food  consumption  for  each  target  species. 

b.  Spatial  analysis  of  fish  biomass.  To  complement  the  biomass-based 
modeling,  we  used  correspondence  analysis  (CA)  (  Legendre  and 
Legendre  1998)  to  examine  regional  variation  in  fish  commimity 
biomass  (standing  crop).  CA  is  an  indirect  gradient  analysis  useful  in 
analyzing  species  x  site  data  matrices.  Eigenvalues  and  axis  loadings 
were  calculated  with  PC-ORD  Version  4  (McCune  and  Mefford  1999). 
Although  CA  has  been  criticized  for  creating  the  spurious  “arch”  effect 
on  the  second  and  subsequent  ordination  axes  (Hill  and  Gauch  1980),  a 
similar  analysis  correcting  for  this  arch  (detrended  CA)  yielded  patterns 
very  similar  to  CA.  Because  rare  species  can  potentially  be  problematic, 
species  that  occurred  only  once  in  &e  data  set  were  omitted  from  the 
analysis;  lumped  taxa  (e.g.,  minnows)  were  also  omitted.  Species 
biomass  estimates  were  square-root  transformed  to  reduce  the 
dominating  effect  of  the  most  common  species  (M.W.  Palmer, 
Ordination  Methods  for  Ecologists  website,  Oklahoma  State  University). 


Model  Sensitivity  Anaiyses 

Sensitivity  analyses  were  used  to  assess  reliability  of  predictions  and  the 
relative  importance  of  model  parameters.  Variability  of  model  output  was 
simulated  with  probabilistic  (Monte  Carlo)  assessment  methodology  wherein 
values  of  model  parameters  were  randomly  selected  fi-om  statistical  distributions. 
Sensitivity  analyses  of  the  food  consumption  model  were  generated  by  altering 
the  five  primary  model  parameters  (water  temperature,  stomach  fullness,  and  the 
intercept  and  slope  coefficients  fi’om  the  gastric  evacuation  model)  by  +  5, 10, 

15, 20,  and  30  percent.  Each  parameter  was  varied  individually,  with  die  others 
held  to  their  original  values.  Random  values  of  each  adjusted  parameter  were 
drawn  firom  wiAin  triangular  constraints  of  each  as  described  earlier.  Similarly, 
sensitivity  of  the  biomass-based  model  involved  varying  the  four  primary 
parameters  (intercept  and  slope  coefficient  of  the  biomass-latitude/longitude 
regression,  P-B  ratio,  and  GCE)  by  +  5, 10, 15, 20,  and  30  percent.  Each 


Chapter  2  Methods  and  Materials 


9 


parameter  was  varied  individually,  with  the  others  constant.  Following  each 
simulation,  the  resulting  percent  change  in  annual  consumption  was  then 
determined  and  used  as  an  index  of  the  model’s  sensitivity  to  that  parameter. 
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3  Results 


Food  Consumption  Model 

Thermal  differences  affecting  rates  of  food  consumption  existed  among  study 
areas  (Figure  2).  In  northern  Lake  Michigan,  water  temperatures  did  not  exceed 
20  “C,  and  were  >15  “C  only  during  a  2-month  period  from  mid- July  through 
mid-September.  The  growing  season  for  warmwater  fishes  (period  with 
temperatures  >  15  °C)  was  estimated  to  be  about  75  days  yeaf  *  (21  percent  of  the 
year)  in  northern  Lake  Michigan  (Table  3).  At  other  end  of  the  spectrum  in  the 
LMR,  water  temperatures  exceeded  20  "C  during  41  percent  of  the  year;  overall 
in  the  LMR,  the  average  growing  season  was  estimated  to  be  200  days  year  . 

The  rank  order  of  warmest  to  coolest  systems  did  not  strictly  correspond  with 
latitudinal  gradient.  The  Ohio  River,  which  flows  approximately  east  to  west  and 
drains  the  Mid- Atlantic  Highlands  region,  was  on  average  warmer  than  the  UMR 
Pool  26,  despite  greater  latitude.  The  Fox  River  and  LTMR  Pool  4  were  also 
warmer  than  Lake  Erie.  In  spite  of  its  latitude.  Lake  Erie  was  the  second  coolest 
system  with  the  second  shortest  growing  season — only  34  percent  of  the  year 
(Table  3).  Lake  Erie  is  the  downstream  outlet  for  all  of  the  upper  Great  Lakes 
(i.e.,  Superior,  Michigan,  and  Huron),  which  are  all  relatively  colder  systems 
supporting  predominately  coolwater  (e.g.,  esocids  and  percids)  and  coldwater 
(e.g.,  salmonid)  fisheries  (Jude  and  Leach  1999). 

T  atitiiHinal  variations  in  atmual  food  consumption  paralleled  water 
temperature.  With  ages  3  to  4  common  carp,  annual  consumption  in  the  LMR 
(warmest  system)  exceeded  that  in  northern  Lake  Michigan  (coolest  system)  by 
117  percent  in  the  average-growth  simulation  (Figure  3).  Similarly,  proportional 
differences  were  1 12  percent  during  the  slow-growth  simulation  and  115  percent 
during  the  fast-growth  simulation  (Figure  3).  Comparable  patterns  were 
observed  with  ages  5  to  6  and  ages  7  to  8  cohorts.  Ranked  from  greatest  to  least, 
annual  food  consumption  ordered  as  follows:  the  LMR,  Ohio  River,  UMR  Pool 
26,  UMR  Pool  13,  UMR  Pool  4,  Fox  River,  Lake  Erie,  and  northern  Lake 
Michigan.  This  ranking  of  annual  food  consumption  followed  exactly  as  did 
mean  annual  water  temperatures  and  nearly  corresponded  with  estimated  len^h 
of  growing  seasons  (Table  3).  Identical  patterns  were  observed  from  simulations 
with  the  other  target  species  (Figures  4  through  6). 
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Figure  2.  Annual  thermographs  for  several  U.S.  river  and  lake  systems  where  zebra  mussels  are 

present.  River  and  lake  systems  are  arranged  A-H  by  increasing  latitude.  Data  provided  by 
the  U.S.  Army  Corps  of  Engineers  and  U.S.  Geological  Survey 
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Table  3 

Thermal  Characteristics  of  the  Eight  River  and  La 
Used  for  Food  Consumption  Simulations.  System 
Latitude 

kes  Syste 
s  Orderei 

System 

Mean  Annual  Water 
Temperature  and  Range 

rc) 

Length  of 

Growina 

Season^ 

%  Days 
S15"C 

%Days 

S20X 

Lower 

Mississippi  River 

16.9  (4.8-29.0) 

200 

55 

41 

Upper  Miss  River 
Pool  26 

13.9  (1.6-27.6) 

167 

46 

34 

Ohio  River 

15.4  (2.4-28.4) 

179 

49 

38 

Upper  Miss  River 
Pool  13 

12.4  (0.0-26.3) 

156 

43 

30 

Lake  Erie 

10.0  (0.0-22.8) 

125 

34 

19 

Fox  River 

11.2  (-0.3-25.8) 

149 

41 

26 

Upper  Miss  River 
Pool  4 

11.4(0.0-27.3) 

160 

44 

28 

Lake  Michigan 

8.2  (2.0-19.3) 

75 

21 

0 

I _ Means 

151 

42 _ 

27 

^  Equivalent  to  the 

number  of  days  during  the  year  with  si  5"C  _ _ _ | 

Size  differences  among  species  directly  contributed  to  observed  differences 
in  food  consumption.  For  instance,  consumption  by  ages  3  to  4  lake  sturgeon 
(greatest  annual  consumption)  exceeded  that  of  ages  3  to  4  channel  catfish  (least 
annual  consumption)  by  2,600  percent  (Table  4).  Similar  trends  in  sturgeon 
consumption  were  observed  over  common  carp  (780  percent)  and  freshwater 
drum  (1,260  percent).  Consumption  by  common  carp  exceeded  that  of 
freshwater  drum  and  channel  catfish  by  55  and  207  percent,  respectively 
(Table  4).  Freshwater  drum  consumption  exceeded  that  of  channel  catfish  by 
99  percent  (Table  4).  Lake  sturgeon  are  large,  long-lived  species  (Billard  and 
Lecointre  2001)  and  mean  size  between  ages  3  and  4  was  estimated  to  be  7.3  kg 
(Carlander  1969).  Thus,  using  1-percent  stomach  fullness  as  a  guideline  in  the 
modeling,  approximately  507  g  of  food  consumption  was  estimated  for  each  day 
in  the  LMR.  These  results  alone  suggest  that  lake  sturgeon  individually  may 
have  the  greatest  potential  to  impact  zebra  mussels  in  the  aquatic  systems 
modeled,  because  of  their  fast  growth  (relative  to  the  other  target  species)  and 
great  size.  However,  other  factors  such  as  the  species’  distribution,  population 
density,  size  structure,  and  local  abiotic  conditions  will  dictate  the  overall 
potential  of  a  specific  population  to  impact  zebra  mussels  in  a  particular  system. 
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tion  (grams  x  1 ,000) 


Low  Growth  DAvg  Growth  H High  Growth  | 

Age  3  to  4 


Age  5  to  6 


100 
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Simulated  annual  food  consumption  of  common  carp  from  eight  U.S.  river  and  lake 
systems  during  a  year  of  average  thermal  conditions.  Simulation  determined  from 
ages  3  to  4,  ages  5  to  6,  and  ages  7  to  8.  River  and  lake  systems  arranged  by  latitude 
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Annual  Food  Conumption  (grams  x  1 ,000) 


Figure  4.  Simulated  annual  food  consumption  of  lake  sturgeon  from  eight  U.S.  river  and  lake 
systems  during  a  year  of  average  thermal  conditions.  Simulation  determined  from 
ages  3  to  4,  ages  5  to  6,  and  ages  7  to  8.  River  and  lake  systems  arranged  by  latitude 
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Simulated  annual  food  consumption  of  freshwater  drum  from  eight  U.S.  river  and  lake 
systems  during  a  year  of  average  thermal  conditions.  Simulation  determined  from 
ages  3  to  4,  ages  5  to  6,  and  ages  7  to  8.  River  and  lake  systems  arranged  by  latitude 
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Figure  6.  Simulated  annual  food  consumption  of  channel  catfish  from  eight  U.S.  river  and  lake 

systems  during  a  year  of  average  thermal  conditions.  Simulation  determined  from 
ages  3  to  4,  ages  5  to  6,  and  ages  7  to  8.  River  and  lake  systems  arranged  by  latitude 
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Table  4 

Mean  Proportional  (%)  Increases^  in  Annual  Consumption  by 

Target  Species  Averaged  Across  Eight  River  and  Lake  Systems. 
Three  values  in  each  cell  represent  means  for  ages  3  to  4,  5  to  6, 
and  7  to  8,  respectively 

Species 

Lake  Sturgeon 

Common  Carp 

Freshwater  Drum 

Common  carp 

780, 1171.951 

Freshwater  drum 
- - — 

1260, 1783. 1184 

54.7,  48.3.  22.2 

I  Channel  catfish 

2600, 2232, 1206 

207,  83.5,  24.3 

98.6,  23.8, 1.7 

I  Percentage  increases  expressed  for  species  listed  in  columns  2  through 
1  listed  in  rows  1  through  3.  For  example,  lake  sturgeon  consumption  excei 

I  consumption  by  780  percent,  1,171  percent,  and  951  percent  for  ages  3 1 

II  respectively. 

4  relative  to  species 
eded  common  carp 

0  4,  5  to  6,  and  7  to  8, 

Interspecific  differences  in  annual  food  consumption  varied  among  cohorts. 
Annual  consumption  by  lake  sturgeon  exceeded  all  other  species,  exceeding 
channel  catfish  by  2,600  percent  for  ages  3  to  4, 2,232  percent  for  ages  5  to  6, 
and  1,206  percent  for  ages  7  to  8  (Table  4).  However,  proportional  differences 
between  the  other  species  were  much  less  because  of  more  similar  sizes  of  the 
cohorts.  Consumption  by  common  carp  exceeded  that  of  freshwater  drum  by 
55  percent  for  ages  3  to  4, 48  percent  for  ages  5  to  6,  and  22  percent  for  ages  7  to 
8  (Table  4).  Similarly,  carp  consumption  exceeded  channel  catfish  by 
207  percent  for  ages  3  to  4,  84  percent  for  ages  5  to  6,  and  24  percent  for  ages  7 
to  8  (Table  4).  Freshwater  drum  consumption  exceeded  that  of  channel  catfish 
by  99  percent  for  ages  3  to  4,  24  percent  for  ages  5  to  6,  and  2  percent  for  ages  7 
to  8  cohom  (Table  4).  Excluding  lake  sturgeon,  differences  in  annual  food 
consumption  among  species  became  less  as  size  increased.  In  other  words,  as 
fish  size  increased,  different  species  exhibited  more  similar  potentials  to  impact 
zebra  mussels.  This  observation  underscores  that  fish  population  size/age 
structure  may  be  an  important  determinant  in  a  species’  ability  to  impact  zebra 
mussels  by  predation. 


Biomass-based  Food  Consumption 

^  The  data  set  used  to  generate  biomass-based  food  consumption  contained  few 
estimates  fi-om  north  of  the  Ohio  River,  Missouri,  and  Kansas  and  many  from  the 
central  and  southern  United  States.  Thus,  biomass-based  results  need 
interpretation  within  that  context.  Total  fish  biomass  ranged  from  14  to 
1,771  kgha-'  and  averaged  333  kgha'  (±13.3)  for  all  systems  in  the  data  set. 
Multiple  linear  regression  modeling  biomass  as  a  function  of  latitude,  longitude, 
and  latitude-longitude  interaction  were  usually  significant,  but  few  models  had  ' 
significant  terms.  Some  significant  latitudinal  and  longitudinal  relationships 
were  detected  with  total  fish  biomass  and  biomass  of  spotted  sucker,  river 
carpsucker,  and  redhorses.  Total  fish  biomass  exhibited  a  slight  inverse 
relationship  with  latitude  (P  =  0.024)  and  direct  relationship  with  longitude 
(P  =  0.047)  (Figure  7),  suggesting  greater  total  fish  biomasses  in  central  U.S. 
systems  (e.g.,  Missouri,  Kansas,  and  Illinois).  Similarly,  biomasses  of  spotted 
sucker  (P  =  0.035  -  0.045)  and  river  carpsucker  (P  =  0.034  -  0.05)  exhibited 
direct  relationships  with  latitude  and  longitude,  which  also  suggested  greater 


18 


Chapter  3  Results 


biomasses  of  those  species  in  central  U.S.  systems  (Figures  8  and  9).  Redhorses 
(P  <  0.0001)  exhibited  significant  relationships  with  latitude  and  longitude 
(P  <  0.0001),  but  the  relationship  was  direct  for  latitude  and  inverse  for  longitude 
(Figure  10).  Results  suggested  greater  biomasses  of  redhorses  in  the  upper  Ohio 
River  drainage  and  Atlantic  seaboard  states  (e.g..  West  Virginia,  Virginia,  and 
North  Carolina)  and  lower  biomasses  in  souAwestem  systems  (e.g.,  Texas, 
Oklahoma,  and  Louisiana).  In  contrast,  despite  that  overall  models  were  highly 
significant,  biomasses  of  common  carp  (P  =  0.47  -  0.54),  freshwater  drum 
(P  =  0.43  -  0.61),  and  channel  catfish  (P  =  0.17  -  0.35)  exhibited  no  significant 
relationships  with  either  latitude  or  longitude  (Figures  1 1  tol3).  Overall 
regression  models  were  not  significant  for  blue  catfish  (P  =  0.19)  or  buffalo 
fishes  (P  =  0.88)  (Figures  14  andl5). 

In  general,  regression  models  indicated  that  the  biomass  of  several  key 
species  likely  to  consume  zebra  mussels  (e.g.,  channel  catfish,  common  carp,  and 
freshwater  drum)  did  not  vary  appreciably  across  latitudes  (10  to  30  kg'ha 
change),  being  similar  in  both  southern  and  northern  systems.  Thus,  these 
species’  ability  to  impact  zebra  mussels  via  predation  varied  little  with  latitude 
(i.e.,  across  systems).  Latitudinal  biomass  increases  of  other  species  likely  to 
consume  zebra  mussels  (e.g.,  spotted  sucker  and  redhorses)  were  detected. 
However,  the  rate  of  biomass  increase  versus  latitude  for  these  species  was  small, 
increasing  on  average  from  2  to  4  kg  ha  *  in  more  southern  latitudes  to  around  5 
to  8  kg  ha  *  in  more  northern  latitudes.  Biomass  increases  detected  with  river 
carpsucker  were  more  substantial,  ranging  from  5  kg-ha  in  lower  latitudes  to 
around  20  kg-ha  '  in  greater  latitudes.  Thus,  we  would  expect  their  impact  to  be 
intermediate  between  common  carp,  freshwater  drum,  channel  catfish,  and 
spotted  sucker-redhorse  groups. 

Predicted  food  consumption  from  biomass-latitude/longitude  regression 
models  provided  mixed  results.  Models  indicated  that  population-level  food 
consumption  decreased  by  25  percent  for  channel  catfish,  49  percent  for 
freshwater  drum,  and  75  percent  for  blue  catfish  between  the  LMR  and  northern 
Lake  Michigan,  the  most  southern  and  northern  systems  modeled  (Table  5).  In 
contrast,  consumption  increased  for  common  carp  by  62  percent  (Table  5)  and 
spotted  sucker  by  27  percent  over  this  same  range  of  systems  (Table  6).  More 
drastic  increases  were  observed  with  river  carpsucker  (901  percent)  and 
redhorses  (412  percent)  (Table  6),  but  as  before,  predicted  biomasses  of  these 
species  were  much  lower  overall.  Overall,  biomass-based  models  indicated  that 
cumulative  food  consumption  of  the  eight  target  species  increased  by  69  percent 
between  the  LMR  and  northern  Lake  Michigan,  presumably  the  result  of  greater 
biomasses  of  some  of  these  species  in  more  northern  systems.  However,  at  the 
same  time,  cumulative  food  consumption  of  the  entire  fish  community  decreased 
by  21  percent  in  northern  Lake  Michigan  compared  to  the  LMR.  Thus,  in  spite 
of  the  proportional  increases  exhibited  by  some  species,  the  total  food  consumed 
by  all  fish  species  declined  one-fifth  with  latitude. 


Chapter  3  Results 


19 


Biomass  (kg-ha-1 )  Biomass  (kg-ha-1 ) 


Figure  7 .  Linear  regression  analysis  of  total  fish  biomass  vs.  latitude  and 

longitude  for  428  reservoir  and  river  systems  in  the  United  States. 
Systems  are  located  from  Texas  east  to  Florida  and  north  to  the  Ohio 
River,  Missouri,  and  Kansas.  P  value  represents  significance  level  of 
the  regression 
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Figure  8.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
spotted  sucker  from  428  reservoir  and  river  systems  in  the  United 
States.  Systems  are  located  from  Texas  east  to  Florida  and  north  to 
the  Ohio  River,  Missouri,  and  Kansas.  P  values  represent  significance 
level  of  the  individual  regressions 
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Figure  9.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
river  carpsucker  from  428  reservoir  and  river  systems  in  the  United 
States.  Systems  are  located  from  Texas  east  to  Florida  and  north  to 
the  Ohio  River,  Missouri,  and  Kansas.  P  values  represent  significance 
level  of  the  individual  regressions 
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Figure  10.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
redhorses  from  428  reservoir  and  river  systems  in  the  United  States. 
Systems  are  located  from  Texas  east  to  Florida  and  north  to  the  Ohio 
River,  Missouri,  and  Kansas.  P  values  represent  significance  level  of 
the  individual  regressions 
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Figure  1 1 .  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
common  carp  from  428  reservoir  and  river  systems  in  the  United 
States.  Systems  are  located  from  Texas  east  to  Florida  and  north  to 
the  Ohio  River,  Missouri,  and  Kansas.  P  values  represent  significance 
level  of  the  individual  regressions 
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Figure  12.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
freshwater  drum  from  428  reservoir  and  river  systems  in  the  United 
States.  Systems  are  located  from  Texas  east  to  Florida  and  north  to 
the  Ohio  River,  Missouri,  and  Kansas.  P  values  represent  significance 
level  of  the  individual  regressions 
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Figure  13.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
channel  catfish  from  428  reservoir  and  river  systems  in  the  United 
States.  Systems  are  located  from  Texas  east  to  Florida  and  north  to 
the  Ohio  River,  Missouri,  and  Kansas.  P  values  represent  significance 
level  of  the  individual  regressions 
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Blue  catfish 


Figure  14.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 

blue  catfish  from  428  reservoir  and  river  systems  in  the  United  States. 
Systems  are  located  from  Texas  east  to  Florida  and  north  to  the  Ohio 
River,  Missouri,  and  Kansas.  P  values  represent  significance  level  of 
the  individual  regressions 
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Figure  15.  Linear  regression  analysis  of  biomass  vs.  latitude  and  longitude  for 
buffalofishes  from  428  reservoir  and  river  systems  in  the  United 
States.  Systems  are  located  from  Texas  east  to  Florida  and  north  to 
the  Ohio  River,  Missouri,  and  Kansas.  P  values  represent  significance 
level  of  the  individual  regressions 
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System 


Lower  Mississippi  River 


Upper  Miss  River  Pool  26 


Ohio  River 


Upper  Miss  River  Pool  13 


Lake  Erie  _ 


Fox  River  _ 


Upper  Miss  River  Poo!  4 


Lake  Michigan 


Table  5 

Predicted  Annual  Food  Consumption  of  Selected  Fish  Species  for  Each  Study  System 
Estimated  with  Biomass-Based  Model _ _=__============^^ 


Estimated  Food  Consumption,  as  kg*ha'^  ±  SE  


Common  carp  |  Freshwater  drum  |  Channel  catfish  |  Biue  catfish 


765  (±31)  518  (±37)  314  (±29) _ ^Q6{±Ze) 

978  (±  32)  453  (±  38) _ 269  (±  30) _ 63  (±  38) 

835  (±31)  196  (±37) _ 216  (±30) _ 41  (±  37) 

1105  (±  32)  426  (±  39) _ 254  (±  30) _ 49  (±  39) 

1010  (±32)  111  (±38) _ 218  (±  30) _ 24  (±  38) 

1194  (±32)  328  (±39) _ 240  (±  30) _ 34  (±  39) 

1245  (±  32)  544  (±  40) _ 246  (±  31) _ 49  (±  39) 

1237  (±  32)  264  (±  40)  _  236  (±  30)  I  27  (±39) 


Note:  Lake  sturgeon  excluded  from  this  evaluation  because  too  few  biomass  estimates  existed  for  a  meaningful  anatygs. 


Tables  ^  ^ 

Predicted  Annual  Food  Consumption  of  Selected  Fish  Species  for  Each  Study 

Estimated  with  Biomass-Based  Model _ _  _ 


Estimated  Food  Consumption,  as  kg*ha  ±  SE _ 


System 


Lower  Mississippi  River 


Upper  Miss  River  Pool  26 


Ohio  River  _ 


Upper  Miss  River  Pool  13 


Lake  Erie 


Fox  River 


Upper  Miss  River  Pool  4 


Lake  Michigan 


Buffalofishes 


726  (±41) 


851  (±44) 


1139  (±  42) 


929  (±  45) 


1880  (±  44) 


1121  (±  45) 


848  (±  46) 


1305  (±  46) 


Spotted  sucker 


252  (±  38) 


400  (±40) 


198  (±  39) 


480  (±  41) 


96  (±  40) 


406  (±42) 


827  (±  42) 


319  (±  42) 


River  carpsucker 


100  (±  39) 


289  (±  41) 


201  (±  40) 


495  (±  42) 


622  (±  41) 


774  (±  43) 


725  (±  43) 


1001  (±43) 


Redhorses 


90  (±  33) 


130  (±  34) 


375  (±  34) 


160  (±  35) 


231  (±  34) 


284  (±35) 


109  (±  35) 


461  (±  35) 


Note:  Lake  sturgeon  excluded  from  this  evaluation  because  too  few  biomass  estimates  existed  for  a  meaningftj[analjgis 


Spatial  Variation  in  Fish  Communities 

Correspondence  analysis  (CA)  explained  35  percent  of  the  variation  in  fish 
communities  on  the  first  three  ordination  axes,  which  suggested  strong  spatial 
differences  in  fish  communities.  The  ordering  of  sites  along  the  first  axis  was 
strongly  associated  with  longitude,  whereas  the  second  and  third  axes  were 
associated  with  latitude.  States  of  the  southern  United  States  tended  to  have 
greater  axis-1  and  axis-2  scores,  whereas  more  northern  states  contained  greater 
axis-1  and  lower  axis-2  scores  (Figure  16).  When  plotting  axis-3  versus  axis-1 
scores,  only  the  latitudinal  gradient  was  evident,  with  axis-3  scores  exhibiting  a 
slight  inverse  relationship  to  latitude  (Figure  17). 
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Figure  16.  Scatterplot  of  axis-1  and  axis-2  scores  from  correspondence  analysis 
of  fish  biomass  estimates  from  U.S.  reservoir  and  river  systems. 
Eigenvalues  were  0.165  for  axis  1  and  0.104  for  axis  2,  which  together 
explained  31.8  percent  of  the  variance  in  biomass.  Axis-1  and  axis-2 
scores  were  averaged  by  state  and  presented  in  lieu  of  showing  all 
individual  data  points.  Arrows  indicate  general  direction  of  gradients. 
Species  listed  at  ends  of  axes  represent  the  rank  order  of  greatest 
axis  loadings 


Several  species  were  identified  as  contributing  to  the  ordering  of  sites  along 
ordination  axes.  Positive  loadings  on  the  first  axis  generally  indicated  that 
southeastern  United  States  systems  (e.g.,  Virginia,  West  Virginia,  North  and 
South  Carolina,  and  Georgia)  contained  greater  biomasses  of  yellow  perch,  white 
catfish,  redbreast  sunfish,  pumpkinseed,  and  chubsuckers.  The  result  with  yellow 
perch  seems  counter-intuitive.  However,  yellow  perch  were  widely  introduced 
into  upland  reservoirs  in  the  southeastern  United  States  in  the  1960s  (Clugston, 
Oliver,  and  Ruelle  1978).  Negative  axis-1  loadings  indicated  greater  biomasses 
of  common  carp,  buffaloes,  river  carpsuckers,  channel  catfish,  blue  catfish,  and 
freshwater  drum  in  southern  (e.g.,  Tennessee,  Alabama,  Mississippi,  Texas, 
Kentucky,  and  Arkansas)  and  central  (e.g.,  Missouri,  Kansas,  Oklahoma,  and 
Illinois)  U.S.  systems.  Axis-2  loadings  mainly  corroborated  axis-1  loadings. 
Positive  axis-2  loadings  indicated  greater  biomasses  of  buffaloes,  spotted  sucker, 
freshwater  drum,  blue  catfish,  chubsuckers,  gars,  bowfins,  and  warmouth  in  more 
southern  U.S.  systems.  Negative  axis-2  loadings  indicated  greater  biomasses  of 
common  carp,  river  carpsucker,  and  walleye  in  central  U.S.  systems.  Overall, 

CA  results  suggested  that  central  and  southern  U.S.  systems  contained  greater 
biomasses  of  several  species  (e.g.,  common  carp,  buffaloes,  freshwater  drum, 
channel  catfish,  spotted  sucker,  and  river  carpsucker)  likely  to  consume  zebra 
mussels. 
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Sensitivity  Anaiyses 

Food  consumption  model 

A  sensitivity  analysis  of  the  parameters  used  in  the  food  consumption  model 
indicated  that  water  temperature  and  slope  of  the  evacuation  rate  model  were 
least  sensitive  to  error,  and  that  stomach  fullness,  the  intercept  term  to  the 
evacuation  rate  model,  and  the  water  temperature  exponent  were  most  sensitive 
to  error.  Increases  and  decreases  in  water  temperature  yielded  differences  up  to 
2  percent  in  annual  consumption  at  a  +30  percent  change  (Table  7).  Similar 
changes  in  the  slope  term  of  the  evacuation  model  yielded  differences  up  to 
5  percent  in  annual  consumption.  Differences  were  greatest  for  the  intercept  term 
of  the  evacuation  rate  model,  which  yielded  a  -59  percent  change  in  aimual 
consumption  when  increased  by  30  percent,  and  +146  percent  change  when 
decreased  by  30  percent  (Table  7).  Similar  changes  in  the  temperature  exponent 
term  yielded  +56  percent  difference  in  annual  consumption  when  the  parameter 
was  increased  by  30  percent  and  a  -34  percent  change  when  the  parameter  was 
decreased  30  percent.  However,  because  water  temperature  daily  means  were 
derived  from  at  least  several  years  of  data  (range  4  to  72  years  per  system, 
average  10  years),  large  errors  in  estimation  were  considered  unlikely. 
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Differences  in  stomach  fullness  were  intermediate  to  these  two  parameters 
(Table  7). 


Table  7 

Proportional  Changes  in  Annual  Consumption  Based  on 
Sensitivity  Analysis  of  Food  Consumption  Modeling 


Change 

% 

Water 

Temperature 

Food  Consumption  Model  Parameters 

Evacuation  Rate 
Model 

Temperature 

Exponent 

Stomach 

Fullness 

Evacuation 
Rate  Model 
Intercept 

Evacuation 
Rate  Model 
Slope 

30 

2.1 

64.0 

-59.4 

5.1 

55.7 

20 

0.8 

40.8 

-45.2 

3.3 

33.9 

10 

0.1 

19.5 

-26.0 

1.7 

15.5 

5 

0.05 

9.5 

-14.0 

0.8 

7.4 

0 

- 

- 

- 

“ 

- 

-5 

-0.1 

‘  -9.1 

16.2 

-0.8 

-6.8 

-10 

-0.3 

-17.8 

35.1 

-1.6 

-13.1 

-20 

-0.7 

-33.6 

82.5 

-3.2 

-24.2 

-30 

-1.2 

-47.7 

147 

-4.8 

-33.6 

Biomass-based  model 


Sensitivity  analysis  of  the  biomass-based  model  indicated  that  all  four 
parameters  were  sensitive  to  error.  In  the  case  of  P-B  ratio  and  GCE,  potential 
error  was  considered  unlikely  because  means  and  ranges  of  these  variables  used 
in  modeling  were  based  on  numerous  reported  studies.  The  intercept  term  in  the 
biomass-latitude  regressions  was  most  sensitive  to  error  (Table  8).  A  30-percent 
increase  in  this  parameter  yielded  a  +15 1  percent  difference  in  annual 
consumption.  Similarly,  a  30-percent  decrease  resulted  in  a  -59  percent  change 
in  consumption.  Changes  in  the  slope  parameter  yielded  smaller  differences  in 
food  consumption  (Table  8). 
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Simulations  indicated  food  consumption  by  individual  fishes  was  inversely 
related  to  latitude  and,  in  effect,  thermal  regimes  of  the  river  and  lake  systems. 
This  observation  was  consistent  with  four  representative  species  likely  to 
consume  zebra  mussels,  namely,  common  carp,  freshwater  drum,  channel  catfish, 
and  lake  sturgeon.  While  rates  of  food  consumption  differed  greatly  among 
species,  annual  food  consumption  in  the  LMR,  the  southernmost  system,  was 
over  100  percent  greater  on  average  than  that  predicted  for  the  northernmost 
system.  Thus,  fishes  in  northern  U.S.  waters  appeared  to  have  less  inherent 
potential  to  impact  zebra  mussel  populations  than  those  in  southern  waters  from 
the  perspective  of  the  rates  of  food  likely  to  be  consumed  during  typical  years. 

Latitudinal  differences  shown  in  this  study  may  actually  be  more  substantial 
in  the  wild  than  indicated  by  modeling.  Average  growth  scenarios  were  derived 
from  both  northern  and  southern  systems.  However,  it  is  likely  that  most 
northern  fish  populations  would  tend  toward  the  slow-growth  scenario  (i.e., 
lower  annual  consumption)  whereas  more  southern  populations  would  fall 
somewhere  between  the  average-  and  fast-growth  scenarios  (i.e.,  greater  annual 
consumption).  In  the  case  of  ages  3  to  4  common  carp,  fast-growth  annual 
consumption  in  the  LMR  exceeded  slow-growth  annual  consumption  in  northern 
Lake  Michigan  by  over  600  percent.  Similar  proportional  differences  were 
observed  with  the  other  species  and  cohorts.  Exceptions  would  likely  exist  in 
that  growth  and  metabolic  demands  of  fishes  in  some  northern  rivers  and  lakes 
would  be  more  similar  to  that  observed  in  many  southern  systems.  However, 
generally  speaking,  it  is  likely  that  food  consumption  by  fishes  in  most  southern 
rivers  and  lakes  would  exceed  most  comparable  northern  waters  by  more  than 
100  percent  during  typical  years.  Although  simulated  food  consumption  by 
individual  fishes  was  consistently  greater  at  lower  latitudes,  the  biomass  of  a 
species  in  a  particular  system  determines  its  potential  to  impact  zebra  mussels. 
Multivariate  analysis  indicated  standing  crop  or  biomass  of  fish  communities 
ch^ged  across  both  latitudinal  and  longitudinal  gradients.  Central  and  southern 
United  States  (i.e.,  Missouri,  Kansas,  Oklahoma,  Kentucky,  Tennessee,  and 
^kansas)  fish  communities  tended  to  have  greater  biomasses  of  common  carp, 
ictalurid  catfishes,  buffalofishes,  carpsuckers,  freshwater  drum,  and  redhorses. 

All  of  these  species  have  been  documented  as  zebra  mussel  predators  (French 
and  Bur  1993;  Tucker  et  al.  1996;  Marsden  1997;  Molloy  et  al.  1997;  Morrison, 
Lynch,  and  Dabrowski  1997),  have  likely  potential  to  consume  zebra  mussels 
based  on  mouthparts  (French  1997),  or  have  been  reported  as  molluscivores.  The 
lack  of  comprehensive  biomass  data  from  extreme  northern  United  States 
systems  with  high  zebra  mussel  densities  (e.g.,  Wisconsin,  Michigan,  New  York, 
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and  the  Great  Lakes)  precluded  a  more  definitive  analysis.  However,  systems  in 
these  regions  generally  support  greater  biomasses  of  percids,  esocids,  and 
salmonids  and  lower  biomasses  of  species  like  common  carp,  freshwater  drum, 
catfishes,  and  catostomid  suckers  (i.e.,  ‘rough’  fishes)  than  comparable  systems 
in  the  central  and  southern  U.S.  (Loftus  and  Regier  1972j  Leach  and  Nepszy 
1976;  Lee  et  al.  1980).  Although  total  fish  biomass  may  not  change  appreciably 
across  latitudes,  the  biomass  contributed  by  species  likely  to  consume  zebra 
mussels  may  change.  Thus,  from  the  composition  of  resident  fish  communities, 
more  northern  communities  appear  to  have  less  potential  to  impact  zebra  mussel 
populations  by  predation  than  their  southern  United  States  counterparts.  In 
addition,  greater  biomasses  of  rough  fishes  in  central  and  southern  United  States 
waters  may  be  helping  to  curb  expansion  of  zebra  mussels  into  the  lower 
Mississippi  River  and  its  tributaries. 

Despite  the  multivariate  findings,  biomass-based  estimates  of  food 
consumption  were  inconsistent  with  the  hypothesis  that  lower  bioi^ses  of  fishes 
that  prey  on  zebra  mussels  in  northern  waters  might  have  implications  over 
control  of  zebra  mussel  densities.  Though  predicted  food  consumption  of  whole 
fish  communities  decreased  by  21  percent  between  the  LMR  (sou&emmost 
system)  and  northern  Lake  Michigan  (northernmost  system),  predicted 
cumulative  consumption  of  the  eight  target  species  increased  by  69  percent  on 
average  between  the  same  systems.  Increases  in  food  consumption  predicted 
from  modeling  for  common  carp  (62  percent),  buffalofishes  (80  percent),  spotted 
sucker  (27  percent),  redhorses  (412  percent),  and  river  carpsucker  (901  percent) 
contributed  most  to  this  observation.  Large  proportional  changes  observed  for 
river  carpsucker  and  redhorses  were  based  on  much  lower  levels  of  biomass  and 
consumption  than  those  predicted  for  common  carp  and  buffalofishes. 

Because  few  biomass  estimates  were  available  for  northern  systems,  we 
believe  that  fish  biomass  and  food  consumption  estimates  may  have  been 
overestimated.  Latitudes  for  northern  Lake  Michigan,  UMR  Pool  4,  Fox  River, 
and  eastern  Lake  Erie  were  the  extreme  end  of  the  range  of  latitudes  used  to 
generate  models.  Thus,  biomass-based  food  consumption  may  have  been 
overestimated  for  these  systems  as  a  result  of  imprecise  estimates  of  biomass. 
Though  exceptions  certainly  exist  for  many  northern  lakes  and  rivers,  greater 
biomasses  of  common  carp,  freshwater  drum,  catfishes,  buffalofishes,  and  other 
catostomids  indicated  from  modeling  are  not  generally  known  from  more 
northern  systems  where  zebra  mussels  are  abundant.  In  10  river  systems  in 
Ohio,  Sanders  et  al.  (1999)  reported  that  freshwater  drum  and  channel  catfish 
were  relatively  minor  components  of  the  fish  community  (<5  percent  in  terms  of 
weight),  especially  in  rivers  draining  into  the  Great  Lakes.  Similarly,  Nelson  and 
Smith  (1981)  categorized  both  channel  catfish  and  freshwater  drum  as  minor 
constituents  of  the  fish  community  in  a  Michigan  river.  Conversely,  common 
carp  abundances  were  highly  variable  in  these  same  systems,  comprising  as  much 
as  70  percent  or  as  little  as  5  percent  of  the  total  fish  biomass  (Nelson  and  Smith 
1981;  Sanders  et  al.  1999).  In  the  Great  Lakes,  where  the  greatest  zebra  mussel 
densities  have  been  reported,  fish  communities  are  dominated  by  salmonids, 
percids,  alewife  Alosa  pseudoharengus,  and  smelt  Osmerus  mordax  (Stewart, 
Kitchell,  and  Crowder  1981;  Jude  and  Leach  1999).  Warmwater  fishes  tend  to 
be  a  minor  part  of  the  fish  community,  though  Lake  Erie  may  be  an  exception 
(Leach  and  Nepszy  1976).  Rivers  and  streams  in  Wisconsin  are  typically 
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coolwater  in  most  areas  and  contain  greater  abundances  of  smallmouth  bass 
Micropterus  dolomieu,  sunfishes  Lepomis  spp.,  salmonids,  and  esocids,  though 
some  systems  in  the  southern  region  of  the  state  tend  to  have  more  warmwater 
assemblages  (Wisconsin  Dept.  Natural  Res.,  unpublished  data).  Though  many  of 
our  target  species  are  resident  to  these  areas,  their  biomasses  generally  do  not 
approach  the  levels  observed  in  more  southern  and  central  U.S.  waters. 
Collectively,  findings  from  other  studies  support  our  contention  of  lower 
biomasses  of  species  likely  to  impact  zebra  mussels  in  more  northern  systems, 
though  modeling  only  partially  supports  this  conclusion  because  of  data 
limitations. 

A  useful  application  based  upon  our  modeling  efforts  would  be  to  assess  how 
much  biomass  is  needed  to  maintain  a  given  level  of  food  consumption  across 
different  systems.  It  is  worth  noting  that  food  consumption  modeling  indicated  a 
100-percent  increase  in  the  rates  of  food  consumed  in  the  LMR  compared  to 
northern  Lake  Michigan  (assuming  average  growth).  Thus,  if  food  consumption 
and  biomass  are  assumed  to  be  related  in  a  1 : 1  fashion,  approximately  twice  the 
fish  biomass  would  be  needed  in  northern  Lake  Michigan  to  maintain  the  same 
level  of  population-level  food  consumption  as  in  the  LMR.  This  conclusion 
assumes  average  growth  in  both  systems,  so  actual  estimates  may  be  greater  or 
lesser.  We  compared  predicted  biomasses  versus  required  biomass  using  this  1 : 1 
biomass-food  consumption  ratio  as  a  guideline.  Because  total  fish  biomass 
declined  slightly  with  latitude,  results  suggested  that  northern  Lake  Michigan 
contained  only  about  one-third  of  the  biomass  required  to  maintain  the  same  level 
of  food  consumption  as  in  the  LMR  (Figure  18).  Proportional  differences  with 
the  other  study  systems  were  smaller  in  magnitude  and  directly  related  to 
latitude.  Predicted  biomass  increases  in  northern  Lake  Michigan  for  common 
carp  (62  percent)  and  spotted  sucker  (27  percent)  did  not  approach  100  percent 
(i.e.,  doubling  of  biomass),  though  buffalofishes  (80  percent)  did  (Figure  19). 
Biomass  increases  predicted  for  river  carpsucker  (901  percent)  and  redhorses 
(408  percent)  far  exceeded  100  percent.  However,  biomass  levels  of  these 
species  were  far  less  than  common  carp  or  buffalofishes  and  would  be  of  less 
consequence  in  terms  of  net  predation  impacts  by  the  whole  fish  community. 

Thus,  several  key  species  of  fishes  may  lack  the  required  biomass  in  more 
northern  latitudes  to  impact  zebra  mussels  the  same  levels  observed  in  southern 
U.S.  waters. 

Although  limitations  of  our  data  sets  and  methodologies  certainly  exist, 
previous  studies  lend  support  to  our  findings.  The  role  of  predators  in  freshwater 
communities  is  complex  and  controversial  (e.g.,  Thorp  1986),  but  previous 
studies  indicate  that  predators  may  significantly  impact  dreissenid  populations. 
Studies  in  European  rivers  and  lakes  (e.g.,  de  Nie  1982;  Daoulas  and  Economids 
1984,  Draulens  1984;  Smit  et  al.  1993)  suggested  that  vertebrate  predators  may 
be  responsible  for  keeping  dreissenid  densities  below  3,000  mussels  m'^ 
compared  to  densities  in  the  Hudson  (approximately  30,000  mussels  m'^)  and 
Ohio  (<40,000  mussels  m'^)  Rivers  (Thorp,  Delong,  and  Casper  1998).  Thorp, 
Delong,  and  Casper  (1998)  also  indicated  that  resident  fishes  impacted  zebra  ’ 
mussels  but  that  impacts  were  greater  among  lower-density  populations.  Given 
the  100-percent  or  greater  increase  in  food  consumption  in  the  LMR  compared  to 
northern  Lake  Michigan  (depending  on  the  growth  scenario  used),  the  lower 
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densities  of  zebra  mussels  in  southern  U.S.  waters  may  be  consistent  with  the 
increased  predation  impacts  of  fish  communities  in  those  waters.  Thus,  fishes  in 
more  southern  U.S.  waters  may  have  a  competitive  edge  in  controlling  or 
restricting  zebra  mussel  populations  because  of  the  bioenergetic  advantages  not 
possible  in  more  northern  rivers  and  lakes. 


(dashed  lines) 


Our  results  were  also  consistent  with  previous  in  situ  studies  despite  differing 
approach  and  differing  scales  of  the  evaluations.  Thorp,  Delong,  and  Casper 
(1998)  reported  fi-om  enclosure  experiments  that  fishes  had  a  stronger  effect  on 
zebra  mussels  in  the  UMR  than  the  Ohio  River.  Furthermore,  some  evidence  of 
size-selective  predation  existed  in  the  UMR,  presumably  because  fishes  were 
able  to  more  easily  select  larger  mussels  firom  less  clumped,  lower-density 
populations.  Though  they  reported  that  zebra  mussels  were  collected  from  the 
guts  of  several  resident  species  (e.g,  smallmouth  buffalo,  channel  and  flathead 
catfishes  Pylodictis  olivaris,  and  river  carpsucker),  they  did  not  discuss  the 
possibility  that  differing  fish  faunas  between  regions  might  have  played  a  role  in 
their  results.  Pitlo  (1987)  reported  that  mean  fish  biomass  (as  kg  •  ha  *)  for 
several  pools  of  the  UMR  (Illinois-Iowa  and  Wisconsin-Minnesota,  n  =  25)  were 
84.7  (+18.9  =  1  SE)  for  common  carp,  18.0  (±5.9)  for  river  carpsucker,  7.2  (±3.7) 
for  spotted  sucker,  24.4  (±9.6)  for  buffalo  fishes,  14.4  (±3.0)  for  channel  catfish, 
and  13.0  (±5.9)  for  fi-eshwater  drum.  Although  biomass  data  were  not  available 
explicitly  for  the  Ohio  River  mainstem,  estimates  from  nearby  Kentucky 
reservoirs  in  the  same  basin  were  24.5  (±4.2)  for  common  carp,  0.3  (±0.2)  for 
river  carpsucker,  3.0  (±0.9)  for  spotted  sucker,  10.0  (±4.2)  for  buffaloes, 

6.0  (±1.0)  for  channel  catfish,  and  8.2  (±4.5)  for  freshwater  drum.  Assuming 
biomass  estimates  (and  resultant  biomass-based  food  consumption)  in  the  Ohio 
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River  mainstem  are  at  least  somewhat  comparable  to  those  of  nearby  reservoirs, 
differences  in  the  composition  of  fish  communities,  particularly  large 
omnivorous  species,  may  partially  explain  the  differential  results  between  the 
two  study  systems  used  by  Thorp,  Delong,  and  Casper  (1998). 


An  interesting  observation  suggested  from  our  modeling  was  the  tremendous 
potential  that  may  exist  with  lake  sturgeon  as  a  biological  control  agent  for  zebra 
mussels.  Historically,  lake  sturgeon  were  found  throughout  the  Mississippi  River 
Valley  and  Great  Lakes  region.  Overfishing  and  habitat  degradation  reduced 
populations  to  extinction  in  many  areas  by  the  year  1900.  However,  though 
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sturgeon  stocks  are  still  declining  worldwide  (Billard  and  Lecointre  2001),  recent 
conservation  efforts  in  the  United  States  have  reestablished  viable  lake  sturgeon 
populations  in  some  areas  of  the  Great  Lakes.  Given  its  large  size,  longevity,  and 
fairly  diverse  feeding  habits,  rapidly  growing  lake  sturgeon  populations  may 
begin  to  significantly  impact  zebra  mussels  in  the  near  future.  Other  sturgeon 
species  inhabiting  the  Mississippi  River  may  also  affect  zebra  mussel 
populations,  but  at  lower  levels  than  suggested  for  lake  sturgeon  because  of  their 
smaller  size. 

Other  factors  unrelated  to  fish  predation  are  also  important  in  regulating 
zebra  mussels  in  southern  U.S.  waters.  Zebra  mussel  invasions  throughout  North 
America  have  occurred  more  rapidly  through  interconnected  waterways  than  in 
isolated  natural  lakes  (Kraft  and  Johnson  2000).  Given  that  most  southern 
reservoirs  are  located  on  large  rivers  and  are  connected  through  the  Mississippi 
River  or  its  tributary  rivers  (e.g.,  Ohio,  Missouri,  Arkansas,  Tennessee, 
Cumberland),  southern  U.S.  waters  are  easily  susceptible  to  zebra  mussel 
invasion.  However,  southern  river  and  lake  systems  are  significantly  warmer  and 
more  turbid  which  may  hamper  zebra  mussel  survival  (Thorp  et  al.  1998;  Allen, 
Thompson,  and  Ramcharan  1999).  Zebra  mussels  thrived  in  temperate  northern 
U.S.  waters  where  summertime  water  temperatures  usually  do  not  exceed  27  to 
28  "C  (McMahon  1996).  Aldridge,  Payne,  and  Miller  (1995)  further  indicated 
from  laboratory  experiments  that  survival  without  extensive  mortality  could 
occur  at  temperatures  up  to  32  °C  .  However,  Allen,  Thompson,  and  Ramcharan 
(1999)  reported  from  field  studies  in  the  LMR  that,  although  annual  growth  was 
comparable  to  northern  regions,  the  pattern  of  growth  was  not,  with  shell  growth 
ceasing  and  tissue  condition  declining  during  July  and  August  when  water 
temperatures  approached  29  to  30  °C.  During  these  midsummer  periods,  mussel 
mortality  increased,  especially  for  individuals  >15  mm.  They  concluded  that 
Slimmer  mortality  was  probable  for  zebra  mussels  in  the  LMR  and  that  the  extent 
would  depend  on  the  duration  of  extreme  water  temperatures,  population  size 
stmcture,  and  condition  during  the  preceding  spring.  Zebra  mussels  are  also  less 
tolerant  to  anoxia  than  most  North  American  bivalves  and  may  be  restricted  to 
well-oxygenated  rivers  and  epilinmetic  zones  of  reservoirs  and  lakes  (Ma^ews 
and  McMahon  1994).  Furthermore,  zebra  mussels  have  been  more  sensitive  to 
drawdowns  than  either  Asiatic  clams  (another  exotic)  or  other  unionids  (Tucker 
et  al.  1997).  Thus,  physical  conditions  and  fish  predation  may  work  in  concert  to 
limit  zebra  mussel  abimdances  and  distribution  in  most  southern  waters. 
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